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Abstract 

A system of a metastable phase with several sorts of heterogeneous centers is consid- 
ered. An analytical theory for the process of decay in such a system has been constructed. 
The free energy of formation of a critical embryo is assumed to be known in the macro- 
scopic approach. At first all asymptotic cases are investigated and then a general inter- 
mediate solution is suggested. All approximate transformations are accomplished with 
the corresponding numerical estimates and analytical justifications. This is the first part 
of the manuscript, the second part follows the first one in the e-print archive. 
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1 Introduction 



The kinetics of condensation on various centers which simultaneously exist in the system will 
be constructed here. The case of condensation of a supersaturated vapor into a state of liquid 
droplets seems to be the simplest case among the first order phase transitions. Traditionally 
this case is considered as a model to introduce some new theoretical constructions in description 
of the first order phase transitions. Ordinary the process of nucleation occurs on dust parti- 
cles or on heterogeneous centers. A case of heterogeneous condensation was a first object in 
regular investigations of the first order phase transitions made by Wilson A first technical 
device constructed by Wilson, i.e. the famous Wilson chamber was based on the effect of het- 
erogeneous condensation (see also [|IJ). A problem of the chamber purification is one of actual 
problems in experimental observations of the homogeneous condensation and the heterogeneous 
way of nucleation is the most usual way of the droplets formation in nature. So, the case of 
heterogeneous condensation which will be investigated is very wide spread in nature. 

The case of condensation is the most well investigated example of the first order phase 
transitions. This leading role of the case of condensation was outlined by creation of the 
classical theory of nucleation by Volmer ||, Becker and Doering ||, Zeldovitch |4]] which gave 
for the first time an expression for the rate of nucleation (i.e. for the rate of an appearance of 
new droplets), which will be an elementary brick in further constructions. The careful analysis 
of the classical theory leads to great number of the publications with various reconsiderations 
of classical expressions. Among them one has to notice the account of the internal degrees of 
freedom in the embryo made by Lothe and Pound |J, another modifications made by Reiss, 
Cohen and Katz ||, Reiss [0, Fisher ||. The application of the density functional theory to 
the first order phase transition made by R.Evans and D.Oxtoby 0, D.Oxtoby and D.Zeng [fL0[j 
allowed to put the microscopic (mesoscopic) models for the condensated substance as the base 
for an expression for the free energy of the critical embryo and for the rate of nucleation. In the 
publications of Oxtoby and Talanquer [fLl|j , Reiss, Tabazadech and Talbot ||12|| , Reiss, Ellerby 



and Weakliem |13|] - ]16| the role of the environment around a droplet is carefully analysed. It 
is equivalent to some specific choice of the statistical ensemble. 

Nevertheless one has to notice that there is no perfect coincidence between concrete theoret- 
ical predictionsS and experimental results. So, this question is still opened in the modern state 
of investigations. But the mentioned expressions for the nucleation rate reproduce a qualita- 
tive behavior of experimental results quite adequately. Moreover, a relative deviation between 
theoretical predictions and experimental results is very smooth function of parameters. 

The mentioned publications allowed to start an investigation of the global kinetics of het- 
erogeneous condensation. The qualitative description of the global evolution during the first 
order phase transition was initiated by Wakeshima ]T7| who considered time lags (characteris- 



tic times) for condensation. The characteristic time of formation of the droplets spectrum was 
necessary to ensure a correct experimental definition of the stationary rate of nucleation. After 



17j the interest for the description of the global evolution was growing continuously. 

The nature and the content of aerosol or heterogeneous particles are so various that different 
approaches to describe the process of condensation are required. Speaking about the aerosol 
particles one supposes that these particles are already the supercritical objects of a liquid 
phase. These objects are growing irreversibly and regularly in time while the objects of a liquid 



1 Certainly, one can not check directly some general recipes as to calculate the objects like statistical sums. 
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phase formed on heterogeneous centers have to overcome the activation barrier. The latter can 
be done fhictuationally and the process of the supercritical embryo appearance resembles the 
homogeneous nucleation, but the height of activation barrier is smaller than in the homogeneous 
case. 

Here the supercritical objects of a liquid phase (may be with a heterogeneous center inside) 
will be called as droplets while the objects of a liquid phase with an arbitrary size and nature 
will be called as "embryos" or as "particles of a new phase". Very often there is no need 
to consider the process of the droplets formation because in the system there is a sufficient 
quiantity of aerosol in the system. This case is extracted by the evident simplification that the 
total number of droplets is already known. Sometimes this case is described as a barrierless 
formation of droplets. The description of this situation one can see in [[18] where this case was 
completely investigated. 

When the activation barrier exists then the total number of droplets is unknown and one 
has to determine a number of droplets and their size spectrum and has to solve a complex non- 
linear problem. To solve this problem numerical calculations were presented in ||19|| , I fZOfl . The 
numerical method of calculations presented in []20| allows to to establish in |19[ some dimen- 
sionless combinations which essentially simplifies the numerical procedure. This simplification 
gives a way to fulfill some rather complex numerical calculations. 

The next step in the development of the global evolution description was the sectional model 



21 1 which simplifies calculations one more time and allows to take into account both nucleation 



and coagulation. In this publication the process of coagulation isn't described because the 
probability of this process is very low and it can be observed only a long time after the end of 
nucleation, i.e. the end of the process of intensive formation of super-critical embryos. 

The main difficulty of the nucleation kinetics is a necessity to take into account the exhaus- 
tion of a vapor phase due to a vapor consumption by droplets. This exhaustion diminishes the 
nucleation rate for new imaginary droplets. Sometimes this influence isn't essential and it is 
shown in |22| that this situation is rather wide spread in laminar flows. In |22[ a theory for this 
case was given and some methods to describe the global evolution were presented. The great 
importance of the problem of the vapor exhaustion around the droplet was stressed by H.Reiss 
in where the stationary profiles around droplets were obtained. 

To take this influence into account one has to solve a non-linear problem. The simplest 
way to analyze this problem is to act in frames of some characteristic values of time and space. 
Namely this way was chosen in |2~3 . 



An effect of a heterogeneous centers exhaustion also diminishes the intensity of the new 
droplets appearance. This effect also has to be taken into account. One has to solve this 
problem in a self consistent way. The theory for the condensation on one type of heterogeneous 



centers can be found in ||25|| , where the approximate self consistent solution of the balance 
equations is given. 

Now one has to specify the situation considered in this publication. One has to mention 
that ordinary in the system there exist several types of heterogeneous centers of different na- 
ture. These centers can have different values of the activation barriers under the given value 
of metastability in the system. Sometimes the difference between centers with neighbour char- 
acteristics is so small that one can speak about the quasi-continuous activity of the band of 
heterogeneous centers. But sometimes the difference in barriers heights between centers with 
neighbor characteristics leads to their different kinetic behavior. 
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To show the variety of situations two examples will be given. In the process of condensation 
on the ions the free energy of the critical embryo depends on the sign of electric charge. As it 



is shown in |26] the free energy F of an embryo with v molecules inside can be presented as 

F = -bv + au 2/3 + ciz/ 1/3 + (ca + c 3 )z/~ 1/3 + c In v + const . (1) 

Here and later all energy-like values are expressed in units of the mean thermal energy; 
a, b, Co, cx, C2, C3 are some parameters. This formula is valid for the embryo in the state of 
internal equilibrium, i.e. near the value v c of a number of molecules in the critical embryo 
(corresponding to the maximum of F) and near the value u e of a number of molecules in the 
equilibrium embryo (corresponding to the minimum of F). All values at v = v c will be marked 
by the lower index c and all values at v — v e will be marked by the lower index e. It is necessary 
to notice that a, b, Cq,Cx, C2 don't depend on the sign of the charge q and a value of C3 is propor- 
tional to q. Since v e < v c then the height of activation barrier AF = F c — F e depends on the 
sign of q. So, in presence of radiation one has two sorts of centers (positive and negative) with 
different heights of activation barriers, i.e. with different activities of heterogeneous centers. 

The spectrum of the sizes of heterogenenous centers as the solid balls with a weak interaction 
between the heterogeneous center and a liquid phase leads to the spectrum of activities of 
heterogeneous centers. In the simple model with the passive nuclei one has to add to the 
surface term the number of the molecules imaginary contained in the volume occupied by the 
nuclei. The free energy of the critical embryo is given by 

F = —bv + a(v + ) 2/3 + const , (2) 

3vi 

where r is the radius of a heterogeneous center, vi is the volume occupied by a molecule of a 
substance in a liquid phase and F e = const. A continuous spectrum of sizes of heterogeneous 
centers initiates a continuous spectrum of the activation barriers heights. 

As the result two different situations can take place - the situation with several types of 
centers and the situation with continuous variation of the properties of centers. Any spectrum 
of activation barrier heights can be split into several continuous parts (may be also like 8- 
functions) which will be considered as the "types" of heterogeneous centers. 

Now one has to specify external conditions. Kinetic description of the condensation pro- 
cess implies external conditions to be known. There are two characteristic types of external 
conditions. The first one is the external conditions of the decay type - until some moment the 
system is in the stable state and then rather rapidly due to the action of external forces the 
initial phase becomes metastable. Then the external action stops. This type of conditions is 
deeply connected with typical conditions in the chamber experiments. Investigations of the 
condensation kinetics [117] were started in this situation. 



The second typical conditions are conditions of the dynamic type. The action of exter- 
nal conditions gradually creates metastability but a rapid growth of a vapor consumption by 
droplets at some moment begin to compensate an action of external conditions, which doesn't 
stop. The rate of a vapor consumption grows so rapidly that the power of metastability falls 
and the process of nucleation stops. This situation is described in for a system with one 
type of heterogeneous centers. 

Why the situation of decay on several types of heterogeneous centers requires a special 
consideration? For the condensation on one type of heterogeneous centers one can see a certain 
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analogy between two mentioned types of external conditions. This is explained by the evident 
notation that in both situations the intensive formation of new droplets occurs near some 
characteristic value of the power of metastabilityi. During the nucleation on different sorts 
of heterogeneous centers in the situation of decay the characteristic value of the power of 
metastability is one and the same for all types of centers and in the situation of dynamic 
external conditions the characteristic value of metastability is determined separately for every 
type of centers. This explains why one has to use different methods for different specific 
situations of heterogeneous condensation on many types of centers. In the current publication 
the situation of decay for heterogeneous centers with arbitrary activities will be analyzed. 

One has to study the period of the intensive formation of the droplets, i.e. the nucleation 
period. The further evolution is described analogously to the condensation on one type of 
centers. 

The following physical assumptions to formulate the model will be used: the thermodynamic 
description of the critical embryo; the random homogeneous space distribution of heterogeneous 
centers; the homogeneous external conditions for temperature and pressure; the high (in com- 
parison with one thermal unit) activation barrieiH; the absence of thermal effects. One can see 
that these assumptions are rather natural. The unit volume is considered. 

According to [27] the condensation equations in the general conditions are analogous to the 
case of the free molecular consumption of the vapor and namely this case will be studied. 

A homogeneous character of a distribution of heterogeneous centers implies that in a char- 
acteristic space region there is a big quantity of heterogeneous centers of every type@. The 
characteristic space size of this region is determined as a mean distance of the diffusion relax- 
ation during the whole nucleation period. For the pure free molecular regime of droplets growth 
the characteristic size is the size of a system. 

The total number of the heterogeneous centers is assumed to be constant in time. This is 
also quite natural because the period of nucleation is relatively short in comparison with other 
characteristic times of the condensation process (except the time of relaxation to the stationary 
state in the near-critical region). 

Since the most interesting characteristics of this process are the numbers of droplets on the 
different types of centers, the accuracy of the theory will be estimated by an error in this value. 



2 Kinetic equations 

The total number of heterogeneous centers will be marked by r] to t i, where % corresponds to 
some sort of the heterogeneous centers. The absence of this index points that the formula is 
valid for an arbitrary sort of heterogeneous centers. The numbers of the free heterogeneous 
centers which aren't occupied by the super critical embryos are marked by the value of rji. 

One can define the neighbor sorts of heterogeneous centers as the centers with neighbor 
values of the activation barriers heights. When the difference in the activation barriers heights 
for the neighbor sorts is many times less than one thermal unit one can speak about the quasi 

2 This value is determined by external conditions. 

3 The theory without the heterogeneous activation barrier is much more simple. 

4 For the quasi continuous spectrum of activities the different type means that the difference of the activation 
barrier height is greater than one thermal unit. 



5 



continuous spectrum of activities of heterogeneous centers. The latter doesn't mean that the 
total variation of the activation barriers heights is less than one thermal unit. 

It is convenient to introduce the characteristic of "activity". One can define "an activity of 
a heterogeneous center" as some parameter w which is proportional to a deviation in the height 
of the activation barrier 

AF(w) = AF \ w=0 -Xw (3) 

with some positive parameter A. The base of decomposition (w = 0) is now an arbitrary value 
and will be chosen later to simplify formulas. 

A total number of heterogeneous centers with a given activity w will be marked by t] tot (w). 
Naturally r] tot {w) is a smooth function of w. Then the neighbor sorts of heterogeneous centers 
have neighbor rates of the droplets formation!. 

A density of the molecules in the equilibrium vapor is marked by n^, a density of the 
molecules in the real vapor is marked by n. A power of the metastability will be described by 
the supersaturation 

Every droplet can be described by the linear size p = u 1 ^ 3 . Due to the free-molecular regime 
of the vapor consumption 

■ w 

where a is the condensation coefficient and r is some characteristic collision time easily obtained 
from the gas kinetic theory. 

The following statements can be analytically proved: 

• The main role in the vapor consumption during nucleation is played by the super critical 
embryos, i.e. by the droplets. 

• The quasi stationary approximation for the nucleation rate is valid during the nucleation 
period. 

Justification of the second statement uses the estimate for the times tf of the establishing of the 
stationary state in the near-critical region which can be found in Q (for heterogeneous case the 
consideration is similar). They have to be many times less than duration of nucleation period. 
There may exist some big times tf. They correspond to big values of u c . It means that AjF 
are huge and these sorts of heterogeneous centers are passive and out of intensive nucleation. 

The characteristic time £* necessary as the base for decompositions will be the time of 
beginning of nucleation. The values at the moment t* will be marked by the lower index *. 
One can choose the zero point of the time axis as t* and introduce the frontal size z according 
to 

z= f (ar^dt' . (5) 
Jo 

Until the coalescence eq. (|5|) ensures the growth of z in time and can be inverted 

t(z)= I'rcT^ . (6) 
Jo C{x) 



3 The sort of the droplets means the sort of the heterogeneous centers. 
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All functions of time become the functions of z and the relative size x = z — p can be introduced. 
During the whole evolution the droplet has one and the same value of the variable x. Con- 
sidering the value t(x) as the moment when the droplet with a given x has been formed (as a 
droplet, i.e. it begins to grow irreversibly) one can consider all functions of time as functions of 
x or z. The variables x and z become equivalent. Hence, one can see that the kinetic equation 
in the supercritical region is reduced to the fact that every droplet keeps the constant value of 
x. To reconstruct the picture of evolution one must establish the dependencies t(z) and ((x). 

The argument oo will mark the total values of the characteristics formed during the whole 
nucleation process. Practically immediately after creation of metatsability the value of the 
supersaturation falls down to 

= C(0) - EiVtot %Ve 1 . (7) 

During the period of the essential formation of the droplets one can assume that the value v e j 
is constant and take it at £ = $*. In initial approximation v e j can be taken at ( = C(0). 

One can analytically prove during the period of nucleation on some given sort for the 
variation of the supersaturation 

I C - C* I < ^ , 

J- i 

where 

dAiF(Q , 

' = * dC lc=$ * ■ (8) 

The same is valid for all sorts of centers when Tj is substituted by miriiTi. 

Let /* i be the amplitude value of the distribution fi of sizes of heterogeneously formed 
droplets measured in units of n^. According to the second statement the value of distribution 
equals to the stationary distribution which is the stationary rate of nucleation divided by the 
droplets rate of growth and by the value of rioo. Then the stationary distribution fi can be 
easily calculated by the following known formulas ||: 

= f = -5^1^, = f c ... (9) 

> Ji i/n« A >• 'h — J( I'll i \VJ 

where W + is the number of the molecules absorbed by the embryo in a unit of time, A e z/ is the 
characteristic width of the equilibrium distribution 

v=(u c +u e )/2 

A ei v= ]T exp(-F(i/)) 

u=l 

and A c jZ/ is the half- width of the near-critical region! 

2 l/2 

A °"= | d*F ,1/2 • 

For the majority of types of heterogeneous centers the following approximations of the 
nucleation rates Ji are valid during the nucleation period 

Ji - Mvtot i, $*) exp(T, (C ~ )-^- . (10) 

Vtot i 



3 Here the steepens descent approximation is used. 
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The validity of these approximations can be easily justified for monotonous interaction between 
the center and the molecules of the substance which is weaker or equal than the function 
reciprocal to a space distance. For the centers with another interaction this approximation has 
to be checked directly. 

The total number of molecules in droplets formed on the sort "i" is marked by n^gi. To 
simplify the formulas 9{ = r)i/r]tat » will be used. 

Using the conservation laws for the heterogeneous centers and for the molecules of the 
substance one can get for g, h 9; L the following equationsB 

g i = f* i [\z-x) 3 ew(-rM^)e i dx = G i (Eg j ,Oi) , (ii) 
6i = exp(-/, i— f Z e X p(-rM^)dx) = SQ» , (12) 

Vtot iJo $* ^ 

where /* , = Ji(r] tot i: $*)r/ ^cm^. These equations form the closed system of the condensation 
equations. This system will be the subject of investigation. 

For the quasi-continuous case one can use the value of activity w as an argument instead 
of the lower index for the sort of heterogeneous centers and get more general approximation 
which covers also the dependence over activity 

(C - $*) 

f((((x),w) = fc($*,w) \ w=0 exp(r — ^— -)exp(w\) , (13) 

where 

r W = . (u) 

The dependence of T on w is rather weak and one can put 

r(iw) = r \ w=0 (15) 

for some essential part of activities spectrum. What is the term "essential part of activities 
spectrum" will be clear later. 

With the help of the conservation laws for heterogeneous centers and for molecules of the 
substance one can get for g, 9 

rz gtot 

g(z, w) = /* / (z- x) 3 exp(-T——)8(x, w)dxr] tot (w) exp(wA) = G w (g tot , 9) , (16) 
Jo < P* 



g tot = J dwg(z,w) , (17) 

rz gtot 

9(z,w) = exp(-/ i(; exp(Au;)n 00 / exp(-T-)dx) = S w (g tot ) , (18) 

Jo 

where /» = f{(®*, w = 0). 



7 Here the first statement is used. 
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The droplets size spectrum can be found as 

f(x, w) = f*exp(\w)exp(-T — - — )6(x,w)r] tot (w) , (19) 

f i = f, i ex P (-r i ^-)e i . (20) 

Since the accuracy of the theory is measured in the terms of the error in the droplets number 
Ni, these values are defined as 

Ni = r ho ti(l-Bi(z))=Q i {B i ) . (21) 
The structure of further constructions will be the following 

• At first the asymptotic cases will be investigated. It will be done for two types of het- 
erogeneous centers. The generalization for many types of centers is evident and trivial. 
Here it isn't assumed that don't depend on i. 

• Then the general intermediate situation will be constructed with the help of approxima- 
tion Ti = const for all sorts of centers. One can prove that this intermediate situation 
and the asymptotic cases cover all possible situations. 

• Then the case of the quasi continuous spectrum of the heterogeneous centers activities 
will be considered. Here one can use some more elegant approaches which will lead to 
more compact results. 



3 Asymptotic solutions 

3.1 Formal generalization of iteration method 

At first one has to show that the case under consideration can not be reduced to trivial gen- 
eralization of known approaches. The formal generalization of the iteration method analogous 



to [28] for the condensation on one given sort of heterogeneous centers leads to the following 



equations: 

9% (i+i) = Gi(z29j (/)) , 9i (i+i) = Si(22gj (i)) , N { = Qi{9i q)) , (22) 

3 3 

9i (o) = , 9i (o) = 1 , (23) 

4 

9% (i) = f*i~r , Oi (i) = exp(-/» i^^z) , (24) 

4 Vtot i 

N t(2) (oo)= Vtotl [l-eM-f^ — (Y,^ 1 )~ 1/AA )} » A = / exp(-a; 4 )rfx « 0.9 . (25) 

Vtot i j 4$„ Jo 

The third iteration can not be calculated in the analytical form. 

Let us analyze A^(2)(oo). Assume that ( is fixed and for some i and j 

f* i ^ /* j ■ 
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Let us decrease n tot j and AjF keeping the constant value of /* j ~ ^ exp(— A,F). It is obvious 
that when 7] tot j is small then the total number of the heterogeneously formed droplets coincides 
with the total quantity of heterogeneous centers and goes to zero when n tot j goes to zero. The 
value of gi at the end of the period of the droplets formation on the heterogeneous centers of 
the sort "j" can be estimated as 

yi — i 

where Axj is the width of the size spectrum (of the size distribution function) for the droplets 
of the sort " f . The value of Axj is restricted from above by the value Axj which is the width 
of the size spectrum without any influence of the droplets of the other sorts and without any 
exhaustion of the heterogeneous centers of this sort. Then the influence of the heterogeneous 
centers of the sort % on the process of the condensation on the centers of the sort j becomes 
negligible in the limit r] tot j — > 0. At the same time ( p5[) shows that in the limit r] tot j — > 0, 
/* i = const the influence of the droplets of the sort i doesn't become small. This leads to the 
big error for Nj (oo) in the second approximation. 

One can not get an analytical expression in the third approximation for Ni in frames of the 
standard iteration method and the second iteration gives the wrong qualitative results. This is 
the main disadvantage of the standard iteration procedure. The reason for this disadvantage 
is the following one. Consider the situation with one type of centers. In the case when the 
interruption of the embryos formation is caused by exhaustion of heterogeneous centers the 
error in the value of gi is compensated by the squeezing force of the operator Si. The analogous 
property is absent for the operator Qi in the situation with several types of centers due to the 
cross influence of the droplets formed on different sorts. This shows that the situation with 
several types of centers can not be effectively described by a formal generalization of the already 
known methods. 



3.2 Characteristic lengths 

The direct generalization of the iteration method fails due to the wrong account of the cross 
influence of the droplets formed on different sorts. Nevertheless it allows to get the spectrum 
of the droplets when the cross influence is excluded. 

On a base of the first iterations in the general procedure one can see that for a separate 
process there are two characteristic lengths. The first one is the length of spectrum in the 
situation when there are no exhaustion of the heterogeneous centers (and no droplets of the 
other sort). One can say that the condensation occurs in the pseudo-homogeneous way. For 
this characteristic value 

^ = (t^) 1/4 ■ (26) 

This length is going from the first iteration for g^. The second length is the length of the 
spectrum the spectrum is formed only by exhaustion of heterogeneous centers. Then the width 
of the spectrum is 

SiX = ^i±- . (27) 
This length is going from the first iteration for 0j. 



10 



Practically the hierarchy between Aj£ and 5iX is ensured by the hierarchy between /* , 
and 77* j. The values of Tj are rather (in comparison with /* $) unsensible to the value of the 
supersaturation. Really 

The value d-Fj c /g?C can De estimated from above by the value in the limit of homogeneous 
condensation dF c hom/ d( when the force of interaction between the heterogeneous center and 
molecules of liquid decreases monotonously in space. Since the energy of the solvatationi 
depends on the supersaturation weaker than F c depends, one can neglect the last term of the 
previous equation and get 

dAF dF r hom / x 

— . (29) 

This dependence is rather weak in comparison with a very sharp dependence of /* j on the 
supersaturation. 

Another important fact is the frontal form of a back side of spectrum© in the pseudo- 
homogeneous situation (when Tj really plays an important role)0. The frontal character can 
be seen from 

fi = fni^{~jS^f. S )^) ■ (30) 

A moderate variation of A^x can be caused only by a very big variation of /* j. 
Instead of Six one can use parameter 

h > = ti (31) 

to simplify formulas. 

Now one can directly analyze all asymptotic cases. 



3.3 The case A\x ~ A2X 
3.3.1 Situation hi < 1, h 2 > 1 

In this situation 

• The process of formation of the droplets on heterogeneous centers of the first sort doesn't 
depend on formation of the droplets of the second sort. 

It can be directly seen from the chain of inequalities 

5%x < Aisc ~ A 2 x < 5 2 x . (32) 
So, the process of formation of the first sort droplets is described by the following equalities: 

9i = hi \\z ~ x) 3 ex V (-T 1 ^-)9 1 dx = G 1 (g 1 , 6 1 ) , (33) 
J 

8 Under barrier character of condensation, i.e. when AiF ^> 1. 
9 The front side has evidently the front character. 

10 The pseudo homogeneous situation can be defined as the situation when the centers of condensation remain 
practically free. 
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Figure 1. 



Solution of the system (\Sal) - (\3o\). The behavior of f = exp(— g)6 as function of x is drawn. 



Vtot 1 ■'O <±>* 



Vtot 

The system ([33T) - ( p4|) can be reduced by rescaling to 

gi = (z — x) 3 exp(— gx{x))6\dx 



t>i = exp( — a / exp(— 1 1 — — )ax) 
JO < 3?* 



(34) 



(35) 



(36) 



with some parameter a = finoo^l^/ivtot iT^ 4 ). Solution of the system (^)-(^) is drawn in 
Fig. 1 for different a. 

The system (^) - fl34j) can be solved by the following iterations 



9i (i+i) = Gx(gi (j),0i (i)) , 0i = ^(pi (j)) , Ni (j) = Qi(0i (»)) . (37) 
When for all values of the arguments 



then 

When for all values of the arguments 
then 



Si(wi) < Si(w 2 ) ■ 
Qi(wi) > Qi(w 2 ) ■ 
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When for all values of the arguments 

Wi < W 2 , V i > v 2 

then 

Gi(vt, w-l) < Gi(v 2 ,w 2 ) ■ 
The initial approximations can be chosen as 

9i (o) = , Q\ (o) = 1 , (38) 

which leads to 

9\ (o) < 9i (2) < ••• < 9i (2i) < ••• < 9\ < ■■■ < 9i (2i+i) < ••• < 9i (3) < 9i (l) , (39) 

vi (o) > vi (2) > ••• > m (2i) > ••• >vi> ■■■ > vi (2i+i) > ••• > vi (3) > m w , 

iVi (o) < N t (2) < ... < N t (2i) < ... < Ni < ... < Ni (2l+l) < ... < N x (3) < Nt (d . (40) 

These estimates allow to prove the convergence of iterations. 
The calculation of iterations gives 



4 

z „ n r 



Pi (1) = /* 1-7- ) #1 (i) = ex P(-/* i ? (41) 

4 7? iot 1 



iVi (2) (oo) = Vtot t (l - «p(-/„ i^(T^)- 1/4 /ri /4 ^)) • (42) 

^toi 1 4$* 

Since 



c/.r 

by the simple numerical calculation of AT 1 ( 3 )(cxd) one can see 



' Ni w - iVx (j) |>0 , (43) 



iVi (2) - iVi 



< 0.015 , 



which means that the second iteration is already rather precise^. 

On the base of the iterations one can get some approximations for the supersaturation0: 

C(i+i) = $*-/* 1 l\z ~ xf exp(-r 1 ^)0 1 {l) dx (44) 
Jo < P* 

and come to the second approximation for £ 

rz ^3 2^2 q z q 

C(2) = $*-/* ij (z- x) 3 exp(-Hx)dx = + + — - — + — exp(-tfz)] , (45) 

where 

# = (46) 
?7tot 1 



n The same procedure can be used for the condensation on one sort of heterogeneous centers. 
12 Without the second sort droplets taken into account. 
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a = Oh = 0.5 



a = 1,6 = 2 



a = 1,6= 0.5 




Figure 2. 



Solution of the system (\5(\) - (51). The behavior of f = exp(— #)# as function of x is drawn. 



This expression can be simplified. The value of supersaturation appears in expression for 
the size spectrum f(x) in the form exp(— Ti(( — $„)/<&*) . After the substitution of C(2) into 
this expression one can see that in the case when ( deviates essentially (i.e. when the exponent 
changes) from all terms except the first two ones can be neglected and 

C(2) = $*-^ 3 ^ • (47) 

For the second sort one can obtain the following system of equations 

92 = /. 2 f\z - x) 3 exp(-r / 2(x) + {r] ' ot l/n °° )x3 )6 2 dx = G 2 (g 2 + ( Vtot 1 /n OB )x 3 , 9 2 ) , (48) 
Jo < P* 

9 2 = exp(-/* 2 / exp(-r 2 )dx) = S 2 {g 2 + {rjtot l/n^x ) . (49) 

Vtot 2 Jo ®* 

This system can be rescaled to 

92= {z — x) 3 exp(— g 2 {x) — bx 3 )9 2 dx , (50) 



o 



9 2 = exp(— a / exp(— g 2 (x) — bx 3 )dx) (51) 



with constants a = U^oo^l /(Vtot ^l^) and b = r] tot 2^1 / (^00$* f¥z ) ■ Solutions of this 
system for different values of a and b are shown in Fig. 2. 
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Having introduced 

A 2 = 92 + (Vtot l/n^z 3 , (52) 



one can rewrite the system (g8|)-(gg) as 

A 2 = G 2 (X 2 , 9 2 ) + ( Vtot ,/n^z 3 = G+(X 2 , 9 2 ) , (53) 

9 2 = S 2 (\ 2 ) . (54) 

The operator G 2 has the same properties as G±, G 2 have. All chains of inequalities remain 
valid with the index "2" instead of "1" and the operator G 2 instead of G 2 . Moreover one can 
see that 

^^AO 1 ^'-^' 1 - ' (55) 

which shows that the worst situation for the iterations convergence is t] tot = 0. Even in this 
situation the second iteration is rather precise. It can be seen from investigation of nucleation 
on the first sort of centers. 

Actually one can avoid here the calculations according to such a complex procedure. The 
term r] tot iz 3 /n 00 ensures the characteristic length 

Di = (-^^)V3 . (56) 

1 2Vtot 1 

Since 

D x > eAiX ~ eA 2 x e ~ (2 3) , (57) 
the condensation on the centers of the second sort occurs in the separate way, the influence of 



the first sort of centers here is negligible and one can use the formulas (p3|)-(p7D, (|41|)-(}43D with 
the index "2" instead of the index "1". 



3.3.2 Situation hi>l, h 2 < 1 

Since ~ A 2 x, one can change the numbers of sorts and reduce this situation to the previous 
one. 



3.3.3 Situation h x > 1, h 2 > 1 

To analyze this situation one has to understand why in the separate condensation on one sort 
already the second iteration gives rather precise results. This property is explained by the 
big power 3 in the subintegral expression for gi. Then the droplets of the big sizes near the 
front side of the spectrum are the main consumers of the vapor. These droplets have small (in 
comparison with Ajx) values of the variable x. The cross influence is rather weak and one can 
use the general iteration procedure which gives 

9i (l) = /* %-r , i = 1,2 , (58) 

9i (i) = exp(-/* i-^^z) , (59) 

Vtot 1 
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N> (2 )(oo) = Vtot ,[1 - «p(-/. ^((Jl)/, i + (^t)f* 2)- 1/4 A)] . (60) 

Since the exhaustion of heterogeneous centers is moderate, the precision of iterations resembles 
the homogeneous case and allows the estimate \N { ( 2 ) — Ni\ < 0.15A^. 

3.3.4 Situation h x < 1, h 2 < 1 

Since 

ftac < Aire ~ A 2 x , (61) 

the droplets of the second sort don't act on the process of formation of the droplets of the first 
sort. Since 

5 2 x < A 2 x ~ Aiar , (62) 

the same is valid for the droplets of the second sort. The system is split into parts corresponding 
to the separate processes of the condensation on different sorts. Eq. (P3"|)-(P7|) can be reproduced 
here. But in the case hi <C 1 for all % the vapor exhaustion can be neglected in comparison with 
the exhaustion of heterogeneous centers and one can get precise explicit results 

Tl 

6Ax) = exp(-/* i — —z) , % = 1,2 , (63) 

Vtot 1 

Tl 

fi(x) = f* »exp(-/, i—^z) , (64) 

Vtot i 

C = $ + _ + Ml)^3 ( 65) 

This expression for £ is obtained by the same procedure as that which led to P7|). 

3.4 Case A x x <C A 2 x 

Due to Aix <C A 2 x the droplets of the second sort don't act on the process of formation of 
the droplets of the first sort. The process^ of formation of the droplets of the first sort can be 



described by the iteration procedure presented by eq. (|37|) 



3.4.1 Situation h x < 1, h 2 > 1 

Due to h\ -C 1 the equation for g\ can be simplified 



9i — f*if (z — x) 3 exp(— Hx)dx ~ ^ tot 1 z 3 . (66) 

JO Hoc 



The value of Q\ is given by (p3|), the value of fi(x) is given by 

For the condensation on the centers of the second sort the equations analogous to (H5| 
are suitable. So, one can get here equations fl52D -( [55f) . But in this situation the inequality (|5 



isn't valid and one must calculate iterations. One can choose as initial approximation \ 2 r \ = 
and get 



A 2 (i) = A 2 Z - + , (67) 
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The case A2X ^> A±x is reversal to this case and can be considered by the simple change of the indexes. 
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2 (2) = exp(-/, rexp(-(-^) 4 - (-^ f)dx) , (68) 

ritot 2 jo Aoo 2 x A h ix 

where 

2* = (^-) 1/4 = A 2 x , A hl x = (^L)V3 = ^ . 

1 2j* 2 1 2%ot 1 

Eq. fl68|) can be rescaled to 

#2 (2) = exp(— a / exp(— x 4 — bx 3 )dx) (69) 
Jo 

with parameters a = /* 2^00 2 x /Vtot 2 and 6 = (A^ 2 x) 3 /(A 1 ^a;) 3 . This dependence is 
shown in Fig. 3. 

The value of 2 x has the sense of the spectrum width when the cross influence and the 
exhaustion of heterogeneous centers are neglected. The value of Ah \x has the sense of the 
spectrum width when the vapor consumption by droplets is neglected. 

One can easily prove that 

4~ I ^2 (o - N 2 {j) |> , — — I N 2 w - N 2 {j) |< 
dx drjtot 1 

for i,j > 2. Then it is easy to show that 

I N 2 - N 2 I 

1 — iW ii < 0.015 
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by the calculation of N 2 (2)(°o) and N 2 (3)(oo) at r] tot i = 0. 

The simple approximation can be obtained if one notices that on the base of (|68|) 



6 2 (2) (oo) - ex P [-/, 2A00 2X^Al + (^) 4 )- 1/4 + f (1 + (^) 3 )- 1/3 )] , (70) 

Vtot 2 * L\h 1% * A>h l x 



where 



B = exp(—x 3 )dy 
Jo 



with the relative error less than 0.035. 

The spectrum of sizes of the droplets formed on the centers of the second sort is 

T 2 /*2^ 4 N , F 2 71 totl 3 ( £ Hoc f z ( ( X N4 t X 



fa = f* 2exp( — — )exp(-- z )exp(-/, 2 / exp(-(- ) -(- ) )dx) . 

$* 4 r] tot 2 Jo 2 x A h ix 

(71) 

By appropriate rescaling it can be reduced to 

fa ~ exp(-z 4 )exp(-(^) 3 )exp(-a jf* exp(-x 4 - (^) 3 )dx) (72) 
with two parameters a — f* 2 n ooAoo 2 x/rj tot 2 and b = xx/Aqo 2 x which is drawn in Fig. 4. 
3.4.2 Situation h x > 1, h 2 < 1 

The description of the process of formation of the droplets of the first sort can't be simplified. 
It has been already given in the previous sections. But the process of formation of the droplets 
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of the second sort is rather simple to describe. The supersaturation is determined by the vapor 
consumption by the droplets of the first sort. Then one has the following expressions 

6 2 = exp(-/, 2 -^~ / exp(-T 2 ^-)dx) , (73) 
Vtot 2 Jo $* 

92 = /* 2 \\z - x) 3 exp(-T 2 ^)6 2 dx . (74) 
Jo <P* 

The value of g 2 during the period of the nucleation on the centers of the second sort can be 
estimated as 

92<^TT , 
J- 2 

which is based on 

5 2 x <C A 2 x . 

So, g 2 is negligible. It is necessary to calculate only 6 2 . To calculate 9 2 one can get into account 
that the value of g\ grows so rapidly that for the value of J Z exp(— T 2 ^-)dx one can show the 
following approximation 

rex P (-r 2 |Mdx«ze(i-^)+ /°°exp(-r 2 |^)rfxe(^-i) , 

Jo $* $* Jo 



or 



where ^ is extracted by the condition 



2i 
r 2 



9i{zb) 

and is the Heavisaid function. The last approximation allows to calculate 6 2 analytically. 
3.4.3 Situation h x > 1, h 2 > 1 

Actually, this situation has been already analyzed in description of the situation A\x <C 
A 2 x; hi > 1, h 2 -C 1 ■ Since Aix <C A 2 x one could not effectively use there the condition h 2 <C 1 
because one could not assume that the inequality h 2 ensures the pure exhaustion of het- 
erogeneous centers without any vapor exhaustion and the consideration made earlier couldn't 
been simplified. As the result the previous consideration covers the situation h\ > 1, h 2 > 1- 



3.4.4 Situation h x < 1, h 2 < 1 

From the first point of view it seems that the situation h\ -C 1, h 2 -C 1 has been already 
described. One has to stress that hj < 1 doesn't allow to state that the nucleation process 
is going at the constant value of the supersaturation. For the first sort nucleation one has 
the previous expressions (|63"D, (|34]). The analogous expressions (^), (0) f° r the second sort 
nucleation can be violated. So, the process of nucleation of the second sort droplets can not be 
described on the base of the initial value of the supersaturation. 
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Figure 5. 
Function 9 2 (2)i x ) f or different a. 



The calculation of g 2 isn't necessary and only the calculation of 6 2 is essential. One has 



e 



2 (2) 



exp(-/* 



Vtot 2 JO 



exp((- 



the final value for 9 2 can be given by 



6» 2 ( 2 )(oo) = exp[-/, 2 A h x x B] . 

Vtot 2 



Eq. ( [75]) can be reduced by the simple rescaling to 

&2 (2) = exp(— a / exp(— x 3 )dx) 
Jo 



(75) 



(76) 



(77) 



with parameter a = ^^Ah \xjr] tot 2 . This behavior is drawn in Fig. 5. 

The mono-disperse approximation for the size spectrum of the first sort droplets is based 
on the evident chain of the inequalities 

Axi ~ 5\x Axx <^ A 2 x . 

When the mono-disperse approximation fails then the vapor consumption isn't essential at all. 

Now all imaginary possible asymptotic situations have been studied. It doesn't mean that 
all of them take place for two given sorts of heterogeneous centers and some given substance. 

The second part of the manuscript which can be found in the cond-mat e-print archive at 
the same publishing date will complete the theory. 
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